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Atom lasers based on rf-outcoupling from a trapped Bose- 
Einstein condensate can be described by a set of generalized, 
coupled Gross-Pitaevskii equations (GPE). 

If not only one but two radio frequencies are used for out- 
coupling, the atoms emerging from the trap have two dif- 
ferent energies and the total wavefunction of the untrapped 
spin-state is a coherent superposition which leads to a pulsed 
atomic beam. 

We present results for such a situation obtained from a 
ID-GPE model for magnetically trapped Rb-87 in the F = 1 
state. The wavefunction of the atomic beam can be approx- 
imated by a sum of two Airy functions. In the limit of weak 
coupling we calculate the intensity analytically. 
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I. INTRODUCTION 

Now, as the creation of trapped Bose-Einstein conden- 
sates |l]-||] has nearly become a standard technique, the 
interest in the field has turned to the manipulation of 
trapped condensates. A major goal is the creation of 
coherent atomic beams emerging from Bose condensates. 
Such a device, now called an atom laser, was first realized 
in the MIT-group y. This experiment created a rather 
small number of strong pulses of atoms. Recently two 
groups were able to produce continuous |5|] or at least 
quasi-continuous pi coherent atomic beams. 

After the initial proposals of atom-laser models (see 
references in m ) the theoretical treatment of atom lasers 
has focused on two major lines: the work in |g |l9[ starts 
with the coupled Gross-Pitaevskii equations (GPE) (and 
generalizations thereof like the Hartree-Fock-Bogoliubov 
theory) and analyzes the corresponding solutions either 
numerically or analytically. 

On the other hand, the studies in po|-p5[ concentrate 
on the characterization of atom lasers using master equa- 
tions analogous to the work on optical lasers and calcu- 
late properties like the linewidth of atom lasers. 

There are two different means of coherent outcoupling, 
i.e getting atoms out of atom traps, which are inves- 
tigated both theoretically and experimentally. A radio 



frequency can be used to flip the spin-state of magneti- 
cally trapped atoms to an untrapped state; this so called 
radio- frequency outcoupling Iq] has been used in [H,p|. 
An output coupler based on Raman transitions was first 
described in |Q and implemented by the NIST-group [||. 
In rf-outcoupling, the atomic beam just falls downwards 
whereas two-photon Raman-outcoupling can be arranged 
in such a way that the outcoupled atoms acquire momen- 
tum and the resulting beam may thus be pointed in a 
specific direction. 

Like in our previous work |19| we concentrate on rf- 
outcoupling describing the setup of the Hansch-group M 
by a ID-model based on coupled GPEs. In |19| we were 
able to characterize the experimental output rate in a sat- 
isfactory way, differences were mainly due to the reduced 
dimensionality of the model. Here, we want to investi- 
gate a slightly different setup that was recently realized 
experimentally: instead of only one rf-field Bloch et al. 
p7t used two radio frequencies for outcoupling. Conse- 
quently, the coherently outcoupled atoms have two dif- 
ferent energies which leads to a pulsed, coherent atomic 
beam. We will analyze the output properties of such a 
device depending on the detunings of the radio frequen- 
cies and their difference. 

The paper has the following structure: In Sec. O we set 
the notation by introducing the basic facts about GPEs 
and rf-atom lasers. In Sec. Ill we first show how station- 
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ary falling waves can be described by appropriate com- 
binations of Airy functions. These are used to model the 
output of the two-mode atom laser and to calculate its 
average output rate. The analytical results of this section 
are then compared with numerical solutions of the cou- 
pled GPEs (Sec. [rv| ). It turns out that the visibihty of 
the pulse pattern of the beam can be related to the cor- 
relation function of the trapped Bose gas (here at T = 0, 
see also |^). In Sec. ^ we summarize our results. 



II. COUPLED GROSS-PITAEVSKII EQUATIONS 
WITH TWO RADIO FREQUENCIES 

The properties of Bose condensates of atomic va- 
pors at zero temperature are very well described by the 
Gross-Pitaevskii equation (GPE) (see |pq| and references 
therein). It is a nonlinear Schrodinger equation for the 
order parameter of the system, namely the macroscopic 
wavefunction of the Bose condensate, and accounts also 
for the coherence properties of the condensate. It may 
thus serve as a starting point for considerations in coher- 



ent atom optics with Bose condensates as it was done in 
1^ : the authors generahzed the GPE to a system of cou- 
pled equations that describe the physics of the trapped 
and untrapped states of a Bose gas coupled via an rf-field 
including the interaction due to collisions. 

In the present article, we modify this approach to a sit- 
uation where the trapped atoms are exposed to two radio 
frequencies. Such a setup has already been analyzed ex- 
perimentally [^ . We focus on a trap with Rb-87 atoms 
in the F = 1 hyperfine-manifold. The three Zeeman sub- 
levels with rn G {—1, 0, 1} are described by macroscopic 
wavefunctions Tpm] after applying the dipole approxima- 
tion the system of generalized GPEs reads 
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We assume that all Zeeman levels interact with the 
same s-wave scattering length ag = llOaBohr so U = 
4TTH^aoN/M. \\^{f,t)\\^ = EmlV'm(^*)P denotes the 
total density in the trap divided by the number TV of 
particles. The interaction between the atoms and the 
radio frequency fields with frequencies Wrf,i and Wrf,2 is 
described by 

Wm^m' — 2hQ(^cos{LUrf,it) -I- cos(a;rf,2i)) 

X {(^m.m' + l + (5m,m'-l) , (2) 

where V, = ^FMBoiu-Brf/lv^^) (ffF = 1/2) is the Rabi 
frequency corresponding to the magnetic field strength 
B,f. 

We now define Wrf — (wrf,i +'-^rf,2)/2 to be the average 
radio frequency and ^rf — ^rf,i — '-^rf,2 as the difference 
between the two driving fields. In order to get rid of the 
high frequencies, we apply the transformation ipm{t) -^ 
e~'™"'^f*Vm(0 and skip all terms oc exp {—2imuirit) (ro- 
tating wave approximation). This turns Eq. (m) into 
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with 



Wm,m' (t) = 2hn COS {5rit/2) {Sm,m' + 1 + Sm,m'-l) ■ (4) 

Equation (H) represents a set of equations in three di- 
mensions that has to be solved numerically. This is com- 
putationally very demanding. In earlier work p,0,na,n9l 



a one-dimensional (ID) approach has proven to be quite 
successful. Usually, the traps used in experiments are 
elongated in the horizontal plane. Like in |1^ we choose 
coordinates where the z-axis points downwards, the long 
axis of the trap is denoted by y, the short horizontal one 
by X. Taking gravity into account, the total effective 
potentials in z-direction are then given by 

Vm,cs{z, t) = -m MujIz^/2 + m/iA 

-Mgz + gi^U{zM\ (5) 

where A = hw-^f—Vos denotes the detuning from the tran- 
sitions at the trap center {VoS = — gFMBohi-Boff). Due to 
gravity the output coupling proceeds mainly in the ver- 
tical direction so it is natural to choose this direction as 
the ID axis for our model. 

To obtain the effective interaction strength g^D in 
Eq. (ra) we require the chemical potential in Thomas- 
Fermi approximation in the ID model (for the trapped 
m = —1 state) to equal the one of the full 3D situation. 
This leads to 
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with az = ^Jh/{Mu!z). 

The ID approximation of Eq. (^) then reads 

ih — 1prn{z,t) = ( ~2M9z2 + '^m,cs{z,t) \ Tpm{z,t) 
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In this work, we concentrate on the weak coupling 
regime of Eq. (R) which is characterized hy Q < ujz 
P, [l2|Jl7| . The condensate is hardly affected in this 
regime, it merely changes its overall size due to the loss 
of atoms by the outcoupling process ||lj] . In [^ we have 
pointed out that the coupling between the trapped and 
untrapped states is concentrated spatially to the cross- 
ing points of the effective potentials in Eq. ^. In the 
present case where two radio frequencies are present, one 
expects two such points of maximal coupling. They are 
located at 



^rf,l,2 = y/2{A±dr{/2)/{hlOz)az 
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corresponding again to the crossing points of the effec- 
tive potentials for either the radio frequency Wrf,! or Wrf,2 
(see Fig. g). In Fig. |l] the densities of the states m — —I 
and (F = 1) are shown for a specific situation. The 
condensate is in the m = —1 state, the plot shows a 
snapshot of the density distributions after 5 ms of out- 
coupling (for further computational details see Sec. |rv| ). 
The distribution for the m = state demonstrates clearly 
that the points z,-{^i^2 sue indeed the points of maximum 
outcoupling; if one considers the two radio frequencies 
separately, these points correspond to the classical turn- 
ing points in the effective potentials in Eq. (0) (see |Q9[ ) . 




FIG. 1. Plot of the densities for m = -1,0 (F = 1) in 
the trap for weak outcoupling. The trapped m = —1 state 
(thick hne, scaled down by a factor of 20) shows some small 
excitations on top of the ground state. The untrapped rn — 
state (thin line) exhibits pulses for z > 23 a^ . The dashed 
lines denote the resonance points z^i^\_2 from Eq. (M), here we 
have uj^n,2 = 2-k x (21000 ± 4500) Hz. 



In the weak coupling region we are considering here, 
the state m = \, which is also untrapped, is hardly pop- 
ulated at all. For this reason, its density distribution is 
not shown in Fig. H and we will not discuss it further. 
Nevertheless, it is included in all numerical calculations. 



III. AIRY FUNCTIONS FOR FALLING ATOMS 
AND OUTPUT FLUX 

A. Stationary, falling waves 

Outside the condensate region, the untrapped m — Q 
state feels only the gravitational potential because the 
mean field contribution in Vo_eff is practically zero there. 
Therefore, the wavefunction for the falling particles is a 
linear combination of Airy functions, which are the solu- 
tions of the Schrodinger equation for a linear potential. 
The usual Airy function Ai(^) does not suffice to describe 
a continuous, stationary beam of falling atoms. In anal- 
ogy to free planar waves, where both cosine and sine are 
used to form a traveling wave oc exp{ikx), we have to use 
a special complex solution of the Schrodinger equation 
for a given energy E 

i;E{z, t) - {Ai{-Uz)) - *Bi(-eB(z))} e-^^'/'^ 

= Al(Cis(z))e-^[®(«-(^»+^*/''l. (9) 

The argument ^siz) simply shifts and rescales the z-axis 
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where 1 is a length scale (see e.g. [Ea). The wavefunc- 
tion tjJE{z, t), though being a solution of the Schrodinger 



equation for a homogenous gravitational field, is not nor- 
malizable due to the exponential growth of Bi(^£;(z)) in 
the classical forbidden region z < —E/{Mg). This is to 
be expected, as it describes a situation where particles 
are flowing steadily to z — > oo without a source of parti- 
cles at some finite z. Nevertheless, we are only interested 
in the classically allowed region where the solution is per- 
fectly analogous to a freely falling wave. 

We will later make use of the asymptotic forms (z — > 
oo) oiM{S,) and 9(0 Hi 
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B. Coherent superposition of two falling w^aves 

The two radio frequencies driving the condensate 
should lead to two different energies in the resulting 
atomic beam (see Fig. g). 





1 SJ 1 li 1 


1 1 

V'-l.eff 

A-5rf/2/ 


a; 


^2 ;\ 






El \ 

1 1 1 i; 1 


\i 1 



10 12 



14 



16 18 20 22 24 26 

FIG. 2. Effective potentials Vm.cff and energies E-i.2 of the 
outcoupled atoms for the two detunings A ± (5rf/2. The two 
vertical lines denote the resonance points 2i.f,i,2. 



The energy of the outcoupled beam for a detuning A 
is given by 
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where fi denotes the chemical potential. If we consider 
the radio frequencies Wrf.i and ll)i{,2, we have to use the 
two detunings Ai_2 = ?i^if,i,2 — Kff to get the energies 
El 2 indicated in Fig. El 

Assuming that the outcoupling process proceeds in a 
fully coherent way, the total wavefunction of the outcou- 
pled beam is given by 



1po{z, t) = Af{%l)Ei {z, t) + stf^E^ {z, t)) 
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where Af is an overall normalization constant and s ac- 
counts for the relative strength of the two contributing 
waves. In experiment, one usually measures the prob- 
ability density |'!/'o(-Z)i)P via some imaging technique. 
If the two energies Ei^2 are not too different from each 
other with respect to the mean energy E = {Ei + E2)/2, 
i.e. hSrf/E <C 1, then Eq. (jl^) can be Taylor-expanded 
around E. Taking the asymptotic forms in Eq. (|ll| ) into 
account, we can write IV-'oC-s^: ^)P a-s 
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Apart from the A^^-dependence, Eq. (|lj) has the typical 
form of an interference term with P being the visibility. 
The maxima of this term describe the falling of pulses 
emerging from the trap. They are located at 
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corresponding to a pulse frequency impulse — '5rf/(27r). 
This is exactly in analogy to "mode-locking" of two 
modes in a laser. 

The phase a oi s — \s\e°' appearing in the above equa- 
tions is responsible for a shift of the interference pattern 
i.e. the pulses. Eq. ( p^ ) describes the output of the atom 
laser only after it has reached a kind of stationary oper- 
ation. We therefore use a merely as a free parameter to 
shift the pulses in time. 



C. Rate of outcoupling 

To determine the normalization factor of the wavefunc- 
tion Af in Eq. ( [l3| ) we use the fact that in ID the flux 
density jz equals directly the rate of atoms going through 
a point, jz can be written as 
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Vz being the velocity of atoms in z-direction. Eq. ( |20| ) is 
obtained by considering only the leading terms in S,e{z) 



for z — > oo. Averaging over time in Eq. (n9^ finally leads 
to 
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In ||l2| , Steck et al. gave a formula for the outcoupling 
rate of a three dimensional trap in the weak coupling 
limit. Analogously, this result can be derived for the ID 
situation we are interested in (see pTj). The rate for 
one resonance point Zics (resulting from a detuning A) is 
given by 
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We now assume that the rates at the two resonance points 
•2rf,i,2 add up to the total outcoupling rate. This rate 
must be equal to jz, so finally with the help of Eq. ( [2l| ) 
we arrive at an expression for the normalization 
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where rj|^^ ^ ^ is obtained from the detunings A ± Sr(/2. 



Inspection of Eq. ( [13D shows that \s\ should reflect the 
relative outcoupling rates at the two resonance points. 
We define it via 
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for the visibility. The numerical results in the next sec- 
tion show that this is indeed the correct expression. If 
the frequency difference S^t is very small, one can also 
use 



which leads to 
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which shows that P is related to the correlation function 
of the trapped Bose gas [g7| . 

IV. NUMERICAL RESULTS 

A. Numerical method and observation of pulses 

To solve the ID GPEs in Eq. (0) numerically we prop- 
agate the wavefunctions with a usual split-operator al- 
gorithm using FFT on a ID grid. As starting condition 



we always use the ground state of the trapped m — —1 
state, which we obtain by imaginary time propagation. 

The coupUng of the Zeeman levels via the rf-field is 
implemented in a very straightforward way: if the right 
hand side of Eq. (R) is written as H^, one can split up the 
Hamiltonian H in an obvious way H = T + V{t) + W{t). 
The kinetic and potential parts T and V{t) are taken care 
of by the FFT-split-operator technique. It is sufficient 
to treat the coupling part W{t) separately in the most 
simple way by keeping only the terms up to first order 
in the propagation time. The propagator from time t to 
t + At then reads 



C/(t, t + At) = Uoit, t + At) 1 



^W{t)At\ 
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where Uq is the split-operator propagator for T + V{t). 

The trap parameters are basically taken from pj, 
the trapping frequency for the m — 1 state is lOz = 
2tt X 127 Hz. We often use harmonic oscillator units, 
the length unit is a^ = 0.95 /im, time is measured in 
1/ujz — 1.3 ms. 

To investigate the output in our atom-laser model 
we have propagated the initial wavefunction for up to 
« 30 ms with different coupling parameters. We kept the 
coupling strength rather low to ensure a clear separation 
of the time scales for Rabi oscillations and for the pulses 
predicted in the previous section. We typically use a 
magnetic field strength of Brf = 0.1 mG that is also used 
in the latest experiments in the Hansch group p7[ . This 
amounts to a Rabi frequency of O = 27r x 70 Hz, which 
is less than the trapping frequency Uz and the pulse fre- 
quency Sit (see below Eq. (111)). All calculations were 
done with A^ = 5 x 10^ atoms in the trap. 

In Fig. (H) we show the density distribution of the out- 
coupled state as a function of both space and time. The 
atomic beam emerging from the condensate clearly ex- 
hibits the pulses predicted in section IIIB. 
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FIG. 3. Plot of the outcoupled density. The atoms in the 
free m = state emerge primarily at the two resonance points 
(cf. Fig. hh and then fall downwards (to the right in the plot). 
One can clearly see the pulse trains emerging from the con- 
densate region that extends to z fa 22 a^. The radio frequency 
detunings in this plot are Ai,2 — 271 x (19.5 ± 3.0) kHz. 



To prove that these pulses really follow the trajectories 
in Eq. ([iq). Fig. ^ shows a con tour plot of the data in 
Fig. ^. As we mentioned in Sec 
a in Eq. (|l^ 
this, the pulses show the correct behavior, they start be- 
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Contour plot of the data in Fig. H. The full lines 



0.001,0.002,0.003. 



are the contour lines for |i/;o(z,t)| a.z 



The dashed lines denote the trajectories of the pulse maxima 
fromEq. (|I|). 



B. Output rates and visibility 

In this section we want to compare the numerical re- 
sults to the analytic predictions of Sec. III. 



All results 

were obtained by keeping the frequency a'rf,2 fixed such 
that the corresponding resonance point is at the center of 
the condensate (2rf,2 = 2;sag = g/^1^ A2 — 27rx 16.5 kHz). 
Ai is then given by Ai = A2 -I- (5rf ; we consider only pos- 
itive rf differences (5rf = . . . 27r x 13 kHz. 

The visibility in Fig. n^is rather high, namely P = 0.90. 
If one increases (5rf , the visibility drops down, as can be 
seen in Fig. @. This is due to the fact that it depends 
mainly on the quotient of the wavefunctions at the reso- 
nance points (see Eqs. ( P4| , P5| , P2|) ), which decreases if one 
pushes Zrf,! further to the boundary of the condensate. 

One can now take Eqs. ( p2| , p5 14) to calculate the vis- 
ibility and the norm of the Airy functions from Sec. Ill 
using the values of |'0-i(2rf.i^| either from the Thomas- 
Fermi (TF) approximation [E8[ or directly from the nu- 
merical calculations. The dashed line in Fig. || is obtained 
in this way using the numerical values, outside the con- 
densate it fits the numerical calculations quite well. The 
Thomas-Fermi result (not shown) is only slightly differ- 
ent. The total output rate rl'* +T],'^ can be calculated 

via Eq. ( p3| ) from TV"^; we compared it to the initial decay 
constant of the condensate occupation and found good 
agreement. 
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FIG. 5. Plot of the density distributions of the atom laser 
(see also Fig. hi). This plot is made after 5ms of operation 
with Sif = 27r X 12000 Hz. The right resonance point is in the 
wing of the condensate so the visibility drops to P — 0.61. 
The dashed line is a plot of the function in Eq. ( |l4| ) with 
the normalization of Eq. (psl). The density of the untrapped 
state goes down to zero around z — AOa^ due to absorbing 
boundary conditions used in the calculations. 



After having shown that the output can be character- 
ized by the wavefunction in Eq. (jlj), we want to compare 
further the semi-analytic results of Sec. HI to the numer- 
ics. For different S^f one can plot the visibilities obtained 
either by fitting the data or by calculating P via (ESh us- 
ing the TF- approximation or the values of |^_i(zrf,i,2)| 



from the numerics. Fig. demonstrates the result: the 
TF-approximation (full line) fits the semi-analytic values 
(triangles) quite well. Like in an experiment, we tried 
to fit the density distribution (n^ to our numerical data 
at a fixed time to get values for P. As the figure shows, 
these values are in a reasonable agreement with the other 
results. 
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FIG. 7. Density distribution of the untrapped state for 
(5rf = 27r X 2000 Hz. The pulses have a substructure origi- 
nating from roughly twice the difference frequency Sr{. 



In Fig. ^ we have not shown values for 6i-{ = 27r x 2 kHz. 
The reason for this is a resonance phenomenon in this 
region that leads to a more complicated pulse structure 
(Fig. 0) and prevents us from fitting (|lj) to the data. 
The origin of this resonance is unclear up to now, there 
are indications however that it is related to the creation 
of small particle-like excitations in the condensate via 
a process that is second order in the output coupling 
strength. 

The small wiggles on top of the condensate density dis- 
tribution in Figs. IT] and pi indicate the presence of such 
excitations also for Jrf away from 27r x 2 kHz. They are 
so small for these frequency differences that the assump- 
tion made above Eq. (|2|) that the outcouphng process 
proceeds as if the two radio frequencies were applied sep- 
arately is really satisfied, i.e. there are hardly any pro- 
cesses of second order in the rf-coupling hQ. 
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S,{/{27t) [kHz] 
FIG. 6. Visibility versus radio frequency difference. The 
full line is obtained from the Thomas-Fermi approximation. 
The triangles are calculated from the numerical values of 
^~i{zr!,i,2) at a specific time. The points denoted by cir- 
cles are obtained by fitting the density distribution ( |l4| ) to 
the numerical data at fixed time (after ~ 7.5 ms of outcou- 
pling). The error for these values is smaller than 5% and was 
conservatively deduced from the fitting procedure. 



V. CONCLUSION 

We have shown in this work that a rf-atom laser driven 
by two radio frequencies — a two-mode atom laser — 
gives rise to a coherently pulsing beam of atoms. We 
have given formulas for the output rate of such a device 
and for the visibility of the pulse pattern of the atomic 
beam. These analytic results are valid for atom lasers 
operated at very low temperature (virtually T — 0) and 
with weak output coupling. 

The main limitation of our results comes from dimen- 
sionality: we use a ID model for our calculations instead 
of a full 3D treatment. In ||l^ we have shown that the 
output rate of a rf-atom laser operated with one fre- 
quency can be calculated by the ID model in qualita- 
tive agreement with experiment. Accordingly, apart from 
quantitative agreement of the output rates our model re- 
produces the experimental features quite well p7[ | . 

We did not address the problem of a pump for a true 
cw-atom laser. Instead, we considered rather small cou- 
plings that ensure a continuous and smooth operation for 



about 20 — 50 ms but with slowly decreasing output rates. 
After this time only very few atoms are left in the atom 
trap. The output of a true cw-atom laser should be very 
similar to that of our situation. 

Coherently pulsing beams of atoms are certainly of in- 
terest for applications in quantum and atom optics. The 
setup described in this paper produces beams that always 
point downwards due to gravity. It should be possible to 
use a similar setup with two different two-photon Raman 
transitions to create a pulsed beam that can be pointed 
in any direction by using the momentum kick from the 
absorption-emission process. 

Finally, if one uses more than two frequencies, one 
might create pulsed beams with shorter pulses while 
keeping their repetition rate fixed in analogy to a mode- 
locked optical laser pa] . 
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